Theoretical study of the optical conductivity of a'-NaV205 
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Using finite temperature diagonalization teciiniques it is sliown tliat tiie quarter-filled t-J-V model 
on a trellis lattice structure provides a quantitative explanation of the highly anisotropic optical 
conductivity of the Q'-NaV205 compound. The combined effects of the short-range Coulomb in- 
teraction and valence fluctuations of V-ions determine the main absorption and the fundamental 
gap. Inter-ladder hopping is necessary for the explanation of the anomalous in-gap absorption and 
generation of spectral weight at high energy. The role of valence fluctuations is explained in terms of 
the domain wall excitations of an anisotropic 2D Ising model in a transverse field close to criticality. 
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The physical properties of low dimensional quantum 
systems received considerable attention from both the- 
oretical and experimental sides because of their uncon- 
ventional spin and charge excitation spectra |l|]. In this 
respect, the vanadate a'-NaV205, initially identified as 
an inorganic spin-Peierls (SP) compound, similar to 
CuGeOa, has attracted great interest as a ladder system 
at quarter filling. The original X-ray structure analysis 
PI was pointing in favor of a non-centrosymmetric struc- 
ture, implying two inequivalent vanadium sites. Hence, 
it was concluded that magnetic spin-1/2 V*"*" ions were 
forming one-dimensional (ID) chains well separated by 
nonmagnetic V^^ chains in the (a, 6)-plane. Based on 
this view the insulating properties appear natural, and 
various theoretical studies were performed using dimer- 
ized Heisenberg chains 0] or spin-phonon models 0, 
which provided satisfactory results, e.g. for the suscepti- 
bility. 

However, recent structure studies P-H] at room tem- 
perature showed that the V-ions are all equivalent with 
valence 4.5, which came as a surprise in view of the in- 
sulating properties of this substance. This puzzle was 
resolved by realizing that the molecular orbital state on 
a rung occupied by one electron is a key element of the 
electronic structure. In combination with the large lo- 
cal repulsion U, which confines the formal charges of V- 
ions to +4 and +5, respectively, the insulating behavior 
and the quasi-lD magnetic properties could be explained 



Surprisingly recent NMR-studies Q revealed the ap- 
pearance of two different V-sites below the structural 
transition at Tc = 34 K g, which clearly indicates that 
the transition is not of spin-Peierls type but rather driven 
by charge ordering. There have been several further ex- 
periments showing anomalous features compared to the 
usual spin-Peierls transition. The ratio between spin gap 
and transition temperature 2A/kBTc ^ 6.5 0,^ is 
much larger than the BCS mean-field value 3.52. The 
changes at Tc of the dielectric constant |0| and the ther- 
mal conductivity flj] differ significantly from the case 



of CuGeOa . Moreover the magnetic field dependence 
ATc oc H^ is only about 20% of what is expected for a 
spin-Peierls transition |15| ]. 

Subsequently several theoretical papers analyzed the 
role of nearest and further neighbor Coulomb interaction, 
proposing zig-zag ||l^,|l3l or linear jlj] charge order in the 
ground state. The nature of the excitation spectra of such 
a charge ordered singlet ground state is an open question 
Ijlql . A test of these concepts could emerge from the very 
distinct experimental optical conductivity data |19| , p0[ . 
These experiments show quite different optical interband 
transitions for a- and 6-polarization, i.e., reflecting the 
anisotropic structure of VO-layers. In addition to the 
larger spectral weight observed for a-polarization, the a- 
spectra show weak low energy in-gap absorption p9[ . 

We study the optical response assuming that the dy- 
namics of electrons is described by an extended t-J-V 
model pl| , p2[ including intra- and inter-rung Coulomb 
repulsion for a 2D system with trellis topology shown in 
Fig. n^. While in a previous numerical study |23] zig-zag 
order was imposed, our results show that the excitation 
spectra are not only determined by the interladder cou- 
pling but also by strong intra-ladder valence fluctuations. 
Our study suggests that a'-NaV205 is close to a quantum 
critical point. The t-J-V model is given by: 



H= J2 *y(cl,.5j-,a + i:f.c.) 
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where cl ^ — cj ^{1 — rii^^a-) are constrained electron cre- 
ation operators, n^ = X^o- ^lo-^*.'^' ^^'^ ^« ^^ ^^^ spin-i 
operator. The sums are over all bonds {i,j) between 
neighbors on the trellis lattice and spin a =t,i- The 
various hopping parameters t.y , intersite Coulomb inter- 
actions Vij , and exchange interactions J^ between vana- 
dium neighbors are defined in Fig. y. The exchange 
term is parametrized as Jy = AtfJU, where U is the 
corresponding on-site Hubbard interaction and we as- 




FIG. 1. Schematic structure of the a'-NaV205 

(a, fe)-planes where circle stands for a vanadium site. The 
black sites define the size of the cluster used in the calcula- 
tion. 

sume that U — 4.0 eV. The calculation of (t{oj) has been 
performed by means of a generalization of the Lanczos 
technique for finite temperature (FTLM) [^. In this 
approach the trace of the thermodynamic expectation 
value is evaluated by Monte-Carlo sampling, while for 
the evaluation of the relevant matrix elements the Lanc- 
zos method is used. Since the system is in an insulating 
state, one can restrict the calculation of 17(0-') to the finite 
frequency response given by the Kubo formula 



1 - e^'^'^ 
am{Lj) = Re 



dTe''^^,U{T)jrn), (2) 



where f3 — l/ksT and jm with m{— a,b) denoting a 
component of the current operator 
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Here _R^ is the m-t]\ component of the vector connecting 
neighboring sites i and j. 

First we consider the optical spectra for a system of 
ladders coupled only by the interaction Vxy, i.e., putting 
txy = 0. In Fig. y we compare results for spectra with 
polarization along a and b direction, respectively, for a 
cluster formed by two ladders whose size is shown in Fig. 
p. Periodic boundary conditions are chosen. We use 
ta = 0.4 and h = 0.2 eV |@,|. For the Coulomb repul- 
sion we consider the symmetric condition Va — Vi, — V 
|l6| with V = 0.8 eV and vary Vxy The a-spectra show a 
strong downward shift when increasing the value of Vxy, 
whereas the 6-spectra hardly change. This behavior re- 
flects the different nature of excited states relevant for 
a- and 6-polarization. Polarization along b connects to 
excitations with two and zero electrons per rung, respec- 
tively. The dynamics of these 'doublons' and 'holons' is 
not much influenced by valence fluctuations, since these 
rungs are fully occupied or empty, and also not by Vxy, 
since the interaction is frustrated (Fig. 1). The non- 
trivial Vxy dependence of the a-spectra will be discussed 
below. Here we only emphasize that the a-spectra are 
sensitive to the ratio 2ta/Vb. Indeed, as one can see in 
Fig. ||, if valence fluctuations are weak (small ta), the 




FIG. 2. Optical conductivity aa{u}) (full line) and cTb{u}) 
(dashed line) for a trellis lattice as shown in Fig.l. The pa- 
rameters used are Va = 14 =0.8, ta =0.4, if,=0.2, txy=0. eV 
at a temperature T=300 K. From a) to c) V^y = 0.,0.8, 1.2 
eV, respectively. Inset: Transverse response {T^ Tq)ui for 
a-polarization at q = (0, 0) for the pseudospin model. 



main weight in the a-spectrum is at energy ^ 2Vb, as 
expected for the case of decoupled ladders (see Fig. |^) , 
but different from the result for larger ta- 

The numerical data in Fig. Ep for txy = contains 
already some of the characteristic features of the optical 
experiments, such as the relative peak positions of the 
spectra with a and b polarization and also their different 
spectral weight. Next we consider the role of the hop- 
ping txy along the zig-zag chains of edge sharing VO5 
pyramids. In Fig. results are shown for ta — 0.4, tf, = 
0.2, txy = 0.15 and Va ^ Vb = 0.8, Vxy = 0.9 eV, where 
the txy value was optimized to have the best agreement 
with the experiments given in the inset p9| |. One im- 
portant effect of flnite txy is the generation of additional 
spectral weight on the high energy side of the spectrum. 
We attribute this to the incoherent motion of excitations 
along the a-direction due to the finite txy. Compared to 
Fig. ||c the a-spectrum shows a prepeak at 0.8 eV and a 
peak at 1.3 eV due to inter-ladder excitations, which falls 
in the range of the shoulder in the experimental spectra. 
There is also a weak low energy absorption below 0.3 
eV, reminiscent of the 'charged magnon' absorption ob- 




FIG. 3. Comparison of normalized optical conductivity 
cra(w)/ta for ta = 0.1 and 0.4 eV which documents the ef- 
fect of valence fluctuations; other parameters as for Fig. Ep. 



served in infrared optics in the range lu ~ 0.01-0.15 eV 
|l9| . We believe that these are charge excitations, which 
are activated by spin-charge and inter-ladder coupling, as 
they disappear in the spinless fermion case or when the 
hopping between the ladders is switched off (Fig. gc). 

For 6-polarization the maximum of the absorption is 
at 1.15 eV, i.e., slightly higher than the main peak in the 
a-spectrum. Again there is high energy spectral weight, 
although less pronounced, which extends up to 2.5 eV. 
There is no evidence of low energy excitations in ab{i^), 
consistent with their absence in experiments. 

Temperature effects are small up to 500 K, however 
there is a decrease of spectral weight in the a-polarization 
spectrum for oj < 0.2 eV and T< 200 K. For the chosen 
parameters, the equal-time density correlation functions 
indicates a tendency towards a zig-zag pattern, although 
the charge distribution itself is uniform. 

To get a deeper understanding of the optical proper- 
ties we turn now to a different representation of the low 
energy degrees of freedom. For large U and Va the rel- 
evant subspace of one electron per rung can be repre- 
sented as eigenstates of spin S and pseudospin T opera- 
tors, where the configurations with T^ = ±i correspond 
to the left/right position of the electron within a rung. 
After projecting out the high energy states we obtain the 
following Hamiltonian for a single ladder [|l7|]: 



Hla: 
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-2t,J2T^- + i2V, + -^)J2T,X 



r+b 
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^ J2iSr ■ Sr+t i)(T+r-+, + H.C.) (4) 

r 

where r indicates a composite site formed by the two 
vanadium ions in the rung. The hopping ta acts like a 
transverse field, while the Coulomb interaction between 
electrons on neighbor rungs r and 7'+6 in the ladder is like 
the zz term in the Ising model. The third term is due to 
the exchange between electrons in different rungs along 




FIG. 4. Optical spectra for a-(full line) and 6-(dashed line) 
polarization; the parameters are the same of Fig. ^3 but for 
finite hopping ia;j;=0.15, and V^y = 0.9 eV at T=300 K. The 
inset shows experimental data by Damascelli et al. [hol. The 
ratio r = la/h of a- and 6-intensities integrated from 0.6 to 
2 eV is Ttheor ~ 1.9 and r^xp ~ 2.2, respectively. 



the b direction. The value of A scales between Va + Vb 
and 2ta depending on which scale of energy is larger. 
The coupling between neighboring ladders yields 



H. 



inter — ladder 
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(f ,±|) and Al - 2Vb + K - V,y. In 



with r2 — ri 

this representation (Ta(j^) is given by the transverse re- 
sponse function (T^^Ty)^ with q ^fo, 0) iit„j = 0. The 
latter is shown in the insets of Fig. g and has been eval- 
uated by means of the FTLM using the same parameters 
as for the t-J-V model and assuming that the spin cor- 
relations can be treated as for valence-bond states, i.e., 
(S,.j • Sra) = — |- As one can see the main peak posi- 
tion and its Vxy dependence are properly described. The 
dominant pole structure is consistent with the excitonic 
nature of the optical absorption. The spectral broaden- 
ing in the complete t-J-V model can be attributed to the 
dynamical character of the spin degrees of freedom. 

The strong Vxy dependence of the peak position in the 
a-spectrum cannot be understood if 2ta/Vb is small, i.e., 
when the ground state has a simple zig-zag charge order 
pattern, corresponding to a classical Neel order of pseu- 
dospins. In this case due to the geometry of the inter- 
ladder bonds the flipping of a pseudospin, i.e., moving an 
electron from left (right) to right (left), does not change 
the number of frustrated bonds with energy ~ Vxy, and 
the excitation energy is 2Vf,. This case is sketched in 
Fig. ^, which shows a single flipped spin in an oth- 
erwise Neel ordered arrangement. By flipping a further 
pseudospin in the neighboring ladder, however, e.g. by 
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FIG. 5. Configuration witli one (a) and two pairs of do- 
main-walls (b). Tlie pseudospin pointing to tlie left (right) 
direction indicate a charge configuration with the electron on 
the left (right) position in the rung, respectively. Black (gray) 
stripes show the bond where electrons of neighboring rungs 
which belong to different ladders feel the Coulomb repulsion 
Vxy before (after) flipping one (a) and two (b) pseudospins. 



an optical excitation, one can reach the configuration of 
Fig. |p with one frustrated bond less compared to the 
starting configuration (Fig. |5p,), which leads to a reduc- 
tion of the excitation energy Vxy. Such processes can 
be neglected in the 2ta/Vb — > limit. However, when 
the system is close to its critical point, i.e., 2ta/Vb ~ 1, 
the ground state is dominated by configurations with a 
high density of domain walls. In this case most optical 
excitation processes occur in the neighborhood of already 
existing domain walls. This implies that near critical! ty 
even though the energy of the ground state is not strongly 
affected by V^y, the excitation spectrum undergoes sig- 
nificant changes. The excitation energy is lowered by the 
repulsion V^y as shown in Fig. 0. 

Finally, we address the low-energy features appearing 
in the a-spectra by considering the single chain Ising 
model in a transverse field, i.e., given by the terms oc ta 
and Vb in Eq. (4). The excitation spectrum of this model 
is described by the domain wall dispersion of the form 
2hcos{q), with h = 2ta/Vb 



A, = Vby/l + h^ + 2hcos{q), with h = 2ta/Vb |^,^ 
There are low-energy charge excitations near q ~ tt con- 
nected to the slow dynamics of the domain walls close 
to the critical point, which can however not be optically 
excited. The optical response for a-polarization is de- 
termined by the excitation of g ~ pairs of domain 
walls whose energy only slightly deviates from the static 
charge gap value 2Vb. Yet this supports the idea that 
the low-energy excitations in Fig. || can be attributed 
to these low-energy charge excitations, namely optically 
activated by spin scattering and by the interladder cou- 
pling. This can be corroborated by a simplified version 
of Hamiltonian (4) where the spins are assumed to be 
disordered and static, i.e., by introducing random vari- 
ables for (Sr • Sr+b)- In contrast to the 'charged magnon' 
scenario |19] the energy scale is determined here by the 
domain wall motion. 

We conclude that the optical response at room- 
temperature can be well described by the quarter-filled 
t-J-V model, assuming that the system is close to criti- 



cality 2ta ~ VJ,. Hence, strong valence fluctuations con- 
nected to domain wall dynamics, which have not been 
considered in previous studies, play a crucial role for the 
understanding of the charge ordering transition and the 
charge pattern of the ground state in the low tempera- 
ture phase. The comparison with the experimental data 
also indicates the importance of both the inter-ladder 
hopping txy and the Coulomb repulsion Vxy Finally, we 
note that the nature of the single particle excitations for 
the present model is a completely open issue. Spectral 
distributions obtained recently by ARPES measurements 
are highly anomalous p^, but have been interpreted by 
the charge-spin separation of a doped Heisenberg chain. 
It is natural to expect that the domain wall excitations 
significantly infiuence the character of the single particle 
electron Green's function. 
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